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CALCULATION OF THE COMPLEX ZEROS OF THE FUNCTION P(z) 
COMPLEMENTARY TO THE INCOMPLETE GAMMA FUNCTION.* 

Bt Phiup Fbankun. 

The gamma function of the complex variable may be written as the 
sum of an integral function and the function 

which has the same poles and the same residues as T{z) itself. The fact 
that P(z) = has exactly foiu" complex roots was first stated, although 
with an insufficient proof, by L. Bourguet.f The method of computing 
these four complex roots used here is based on an existence proof given 
recently by T. H. Gronwall.$ This proof consists in first showing that the 
auxiliary rational function 

has all its real roots lying one in each of the 2n — 3 intervals : 

a 1 

— 2m — 1 — TK— vi = 3S— 2m — 1 — ,r, \, (w = 2, 3, • • •, n) 
{2m) I {2m) I ^ > > > / 

-2m-2 + (2m+l)! =' = ~^^~^ + (2m+l)! (^ = 2,3, • • •,n-l), 

and hence that Pn{z) has four complex roots and from this proceeding to 
P{z) itself by appljdng Rouch^'s theorem. The method of computation 
consists in locating the complex zeros of Pio{z) to three decimal places and 
then, using these as initial values, finding the roots of P{z) = by a direct 
application of Newton's method of approximation. 

Starting with the mean values of the above-mentioned intervals for 
n = 10, and applying Newton's method of approximation to 

(1) Pioiz)=t^-^-^ 

* Presented at the meeting of the American Mathematical Society, at New York City, April 
26, 1919. 

t L. Bourguet, Sur la tWorie des int^grales Euleriennes, Comptes rendus, vol. 96 (1883), 
pp. 1487-1490. 

t T. H. Gronwall, Svir les z&os des fonctions P{z) et Q(z) associ^es h, la fonction gamma, 
Annales de I'Ecole Normale, Ser. Ill, vol. 33 (1916), pp. 381-393. 
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we find that the real roots are 

P6 = - 5.11907; ps= - 5.97158; p? = - 7.00432; 

pg = - 7.99940; ps = - 9.00008; pw = - 9.99999; 

p„. m (m = 11, 12, • • • 21), 

where, as in the rest of this computation, quantities which do not affect 
the fifth decimal place are neglected. Hence Pio{z) which may be written as 

Oqs'^ + Oig^ + • • • + OzqS + an 
z{z + l){z + 2) ... (2 + 21) 
is equal to 



(2) 



bog^" + feiz" + • . • + bgg + &10 
z{z + l){z + 2) ...(«+ 10) 



to the desired degree of accuracy. 

If Pio{z) be expanded in powers of 1/z, the coefiicient of l/g' is 



and since this is sensibly equal to 

^(- l)''+'»''' 



v=l Vl 



(r < 10), 



(3) -e{\+h-h^ ?+'■') 

is obtained as a second approximate form for Pio(«). A comparison of (3) 
with (2) gives &o = 1/e, &i = 56/e, 62 = 1375/e. 
The expansion of Pv,{z) in the form 

(4) l + At + A2Z + A^^ + '•• 

(where Ai = SS=i[(- l)''/»'!»'] = .7965996) when compared with (2) deterr 
mines the values h = 7,737,939, few = 3,628,800. 

The numerator of (2) (with the factor 1/e removed) is therefore equal to 

(5) gi" + 562" + 1375z8 + 1- 21,033,9002 + 9,864,102. 

This expression, set equal tjo zero, has for roots the foiu* complex roots of 
Pio{z) = 0, together with ps, pe, • . • pio found above. The product 
IU?=6(2 — Pm) is found to be 

(6) z" + 45.094425 + 355.917«« + • • • + 129,7492 + 154,152. 
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Dividing (5) by (6), we obtain the equation 

(7) 2* + 10.906«8 + 44.4252^ + 82.5892 + 63.989 = 0. 

The determination of the fourth and fifth terms of (7) from the last two 
terms of (5) and (6) is preferable to the use of the first five terms of (5) 
and (6), as it involves fewer operations and consequently gives these terms 
more accurately than the alternative method. 

The roots of (7), obtained by the use of Ferrari's formulae, are, to 
three decimal places: 

Pi = - 1.726 + 1.238i, Pi = - 1.726 - 1.2BSi, 

Ps = - 3.727 + .5441, pi = - 3.727 - .54Ai. 

Starting with these values, and applying Newton's method of approxi- 
mation to P(z) itself, we find that the four complex roots of P(z) are 

2i = - 1.7262976 + 1.2380921i, z^ = - 1.7262976 - 1.2380921i, 

•23 = - 3.7264730 + .5406746i, Zt = - 3.7264730 - .5406746^. 
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